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1 Abstract 

This research is focused on assessing the value of the 
Reynolds- Averaged Navier-Stokes (RANS) method- 
ology for active flow control applications. An exper- 
imental flow control database exists for a TAU0015 
airfoil, which is a modification of a NACA0015 air- 
foil. The airfoil has discontinuities at the leading 
edge due to the implementation of a fluidic actua- 
tor and aft of mid chord on the upper surface. This 
paper documents two- and three-dimensional com- 
putational results for the baseline wing configura- 
tion (no control) with the experimental results. The 
two-dimensional results suggest that the mid-chord 
discontinuity does not effect the aerodynamics of the 
wing and can be ignored for more efficient computa- 
tions. The leading-edge discontinuity significantly 
affects the lift and drag; hence, the integrity of 
the leading-edge notch discontinuity must be main- 
tained in the computations to achieve a good match 
with the experimental data. The three-dimensional 
integrated performance results are in good agree- 
ment with the experiments inspite of some conver- 
gence and grid resolution issues. 

2 Introduction 

Increased attention has been devoted to the develop- 
ment of techniques 1,2 capable of enhancing our abil- 
ity to control steady and unsteady flows in a wide 
variety of configurations and applications. Control- 
ling the flow in and around these configurations can 
lead to improved efficiency and performance com- 
pared to uncontrolled configurations. While pas- 
sive devices, such as the micro vortex generator 3 * 4 
have successfully been incorporated into produc- 
tion vehicles, passive devices are primarily limited 
to a single function or are effective over small op- 
erational ranges. Although decades of research 


have been conducted using classical active controls 
(e.g., suction laminar flow control 5 ), questions of 
system complexity and/or cost versus benefit have 
prevented such controls from reaching marketable 
applications. The recent invention of novel actu- 
ator concepts, some of which are locally applied 
and require only electrical power input, has intro- 
duced the potential for notably new vehicle con- 
trol. Or, these new' actuator /effector technologies 
may cause a paradigm shift in the design of con- 
ventional systems. Such a paradigm shift in high- 
lift system design for transport aircraft has led to 
a projected decrease in the weight, cost, and part- 
card count of a conventional vehicle with flow* con- 
trol augmented high-lift systems. 6 Such benefits are 
then compounded into fuel savings and emission and 
noise reductions. 

A limited number of active flow control appli- 
cations have been tested in the laboratory and in 
flight . These applications include dynamic stall con- 
trol using a deformable leading edge, 7 separation 
control for takeoff and landing flight conditions us- 
ing piezo devices, 8 * 9 pulsed vortex generators, 10 and 
zero-net-mass oscillators, 11-13 and thrust vectoring 
with piezo-driven oscillatory actuation. 11 

Laboratory and flight demonstrations of technol- 
ogy are essential to assess the technical feasibility of 
an innovative concept; however, accurate design and 
analysis tools can potentially lead to revolutionary 
new' configurations and are a critical component in 
the estimation of the cost versus benefit of the tech- 
nology on a real application. Hence, this paper will 
focus on the validation of design tools tow ard active 
flow' control applications. 

As yet, there is no definitive comparison w'ith 
experimental data that indicates current compu- 
tational capabilities can quantitatively predict the 
large aerodynamic performance gains achieved with 
active flow* control in the laboratory. How'ever, one 
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study 15 using the RANS methodology has shown 
good quantitative agreement with experimental re- 
sults for an isolated zero-net-mass actuator. In ad- 
dition, some recent studies 16-18 have used RANS 
to demonstrate qualitative performance gains com- 
pared with the experimental data for separation 
control on an airfoil. Quantitative comparisons 
for both baseline and flow control cases indicated 
that computational results were in poor quantita- 
tive agreement with the experiments. 16,17 

The current research will investigate the potential 
use of an unstructured grid RANS approach to pre- 
dict aerodynamic performance for active flow con- 
trol applications building on the early studies. 16,17 
Before proceeding to the time-dependent flow con- 
trol case, the computational results must quan- 
titatively match experiments for the no control 
case. This paper compares the two- and three- 
dimensional simulation results with experiments for 
the baseline (no control) case. 

The next section describes the configurations used 
for the computations and the experiments. The 
computational approach is then described followed 
by two- and three-dimensional results and conclud- 
ing remarks. 



3 Airfoil Configurations and 
Experimental Data 

Figure 1 shows the leading edge and trailing edge 
regions of the configurations used in the present 
study. The mid-chord regions for all airfoils merge 
to the same configuration at 20% chord and there- 
fore are not shown in figure 1. A NACA0015 airfoil 
is used in the study to provide reference point re- 
sults. The blunt trailing edge of the NACA0015 
airfoil matches the TAU0015. The TAU0015 air- 
foil model was tested in a low-speed wind tunnel 
at the Tel- Aviv University (TAU). The model is a 
NACA0015 airfoil modified in the leading-edge re- 
gion to acommodate an actuation slot. Hence, the 
airfoil configuration is referred to as the TAU0013 
airfoil. The TAU0Oi5 model has a 0.3645 m chord, 
a 0.3% thick blunt trailing edge, and a 0.4% chord 
thick notch at 76.6% chord, which results from the 
flap/main element connection used in a different ex- 
periment. The actuator for the TAU0015 tests was 
located at the leading edge and leads to the 0.3% 
chord discontinuity (straight horizontal line region). 

The previous two-dimensional computational 
studies 16,17 used single-block structured grid RANS 



Figure 1: TAU0015, TAU0015m, and NACA0015 
airfoils. 
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and for convenience ignored tlie 76.6% chord 
notch and trailing edge thickness and faired over 
(smoothed) the leading edge discontinuity. Here, 
this altered configuration is shown in figure 1 and 
is referred to as TAU0015m. Computational re- 
sults for the TAUOOlSm will be compared with the 
TAU0015 configuration. There is no experimental 
data for the TAUOOlSm model. 

The flow control experiments with the TAU0015 
airfoil were conducted in the Meadow- Knapp Low 
Speed Wind Tunnel at the Tel Aviv University. 11 
The test section is 1.50 m high and 0.61 m wide. The 
TAU0015 model was instrumented with 36 static 
pressure taps located on the model centerline and 
measurements were made sequentially using a multi- 
port, mechanical rotary valve (Scanivalve) and a 5 
psi full scale differential pressure transducer. The 
transducer has an accuracy of 0.06% full scale. The 
free stream velocity of all tests was nominally 51 
m/s. Tin 1 pressure coefficient results are accurate to 
within ±0.6%. Lift (C/) and drag (Cd) coefficients 
are obtained by integrating the measured pressures; 
accuracy in Ci is estimated to be ±0.01 for pre- 
stall conditions and ±0.03 post-stall. The drag co- 
efficient Cd has experimental uncertainty of ±0.003 
at pre-stall conditions and and ±10% at post-stall 
conditions. 19 The experimental conditions are at a 
Mach number of 0.15 and a chord Reynolds numbers 
of 1.2 million. The uncertainty in Reynolds number 
is ±2% due to variations in temperature and veloc- 
ity during the tests. 

The Ci and Cd for the experimental data is avail- 
able at angle of attacks (a) from 0° to 24° in 2° in- 
crements. The maximum lift coefficient is Cj = 1.056 
at q = 12°. In addition, pressure coefficients (C p ) 
at a = 8°, a = 12°, a = 14°, and a = 22° are used 
for the current comparison. 

4 Computational Approach 

The fully unstructured Navier Stokes two di- 
mensional (FUN2D) 20 and three dimensional 
(FUN3D) 21 codes solve the time dependent 
Reynolds- Averaged Navier-Stokes equations. The 
equations are expressed as a system of conserva- 
tion laws relating the rate of change of mass, mo- 
mentum, and energy in a control volume to the 
fluxes of these quantities through the control vol- 
ume. The solver is an implicit, upwind-differencing 
algorithm with the inviscid fluxes obtained on each 
face of the control volume, employing Roe’s flux- 


difference-splitting scheme. The node-based algo- 
rithm stores the variables at the vertices of the mesh 
and the equations are solved on the non-overlapping 
control volumes surrounding each node. The vis- 
cous terms are discretized using a central-difference 
formulation. Time-advancement is made with a 
linearized backward Euler scheme. At each time 
step, the equations are solved with 15 Gauss-Seidel 
sub-iterations, sequentially solving for all odd num- 
bered nodes and then all even numbered nodes. 
A two-level V-type multigrid approach is used for 
the two-dimensional baseline study to accelerate 
convergence. 21 No multigrid acceleration was used 
for the three-dimensional calculations. A very fine 
grid was first used to determine the affect of the 
geometric discontinuities on the aerodynamic per- 
formance. Grid coarsening was implemented to im- 
prove computational efficiency without sacraficing 
accuracy. The Spalart-Allmaras turbulence model 22 
is used in this investigation and all computations 
assume fully turbulent flow. A single- processor SGI 
OCTANE was used for the two-dimensional com- 
putations. The three-dimensional calculations were 
conducted on a SGI ORIGIN 2000 machine using 8 
processors. 

The two-dimensional unstructured grids were 
generated with advancing front type point place- 
ment with iterative local re-meshing for grid quality 
improvement. 23,24 Figure 2 shows an initial coarse 
grid around the TAU0015 airfoil, resolving all dis- 
continuities on the configuration. Similar grids were 
generated for the TAU0015m and NACA0015 air- 
foils. The grids extend from the airfoils to form 
a far-field circle with a radius of 20 chord lengths 
around the airfoils. 

The three-dimensional simulations were con- 
ducted on the TAU0015 to validate the computa- 
tional approach with the experiments. Future ac- 
tive flow control simulations will involve the use of 
fully three-dimensional actuation technologies and 
therefore a validation of the three-dimensional code 
on the baseline (no control) cases is essential. 

The three-dimensional TAU0015 grids were 
obtained using the TetrUSS unstructured grid 
package. 25 27 The grid generation technique is 
based on the Advancing-Front (AFM) and the 
Advancing-Layers (ALM) methods and produce 
fully tetrahedral meshes. In defining the compu- 
tational domain, boundaries are represented by bi- 
linear surface patches. The patches are defined 
only in terms of their edges, which are constructed 
in GRIDTOOL based on user- specified geometries. 
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Figure 2: Two-dimensional grid for TAU0015 airfoil 
(top), leading edge region (middle), and aft portion 
of wing (bottom). 


Grid characteristics such as cell spacing and stretch- 
ing are also specified in GRIDTOOL by the place- 
ment of cell "sources”. The surface mesh is gener- 
ated in VGRID by triangulating each surface patch 
with a two-dimensional version of the advancing- 
front method. The triangulated surface patches 
then form the initial "front” for the advancement 
of the three-dimensional tetrahedral volume cells. 
The high aspect ratio cells in the viscous region are 
constructed with the ALM. where individual cells 
are formed by inserting a point based on surface 
vector information. Finally, the isotropic cells were 
constructed in the inviscid region using the AFM. 
Smooth propagation of grid spacing is achieved by 
solving a Poisson equation on a cartesian back- 
ground grid, using the GRIDTOOL defined cell 
sources as inputs. Once the volume grid is gener- 
ated, local remeshing can be implemented for grid 
quality improvements. 

The three-dimensional TAU0015 surface grid has 
a spanwise length of approximately one chord 
length. The outer boundary of the grid extends 10 
chord lengths out from the airfoil geometry form- 
ing a box. Figure 3 shows the surface grid at the 
leading edge of the TALTI015 airfoil. The cells are 
clustered near the leading edge step region to resolve 
the geometry and capture the flow field characteris- 
tics at this discontinuity. The volume grid has 1.85 
million cells. The sidew r alls of the computational 
domain were treated as inviscid walls, with the re- 
maining outer boundaries treated as characteristic 
inflow/outflow surfaces. No- slip viscous boundary 
conditions were prescribed on the airfoil surface. A 
more thorough discussion of the three-dimensional 
grids will be provided in the three-dimensional re- 
sults section. 



Figure 3: Three-dimensional surface grid for 

TAL-0015 airfoil. 
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Case 

Gridl 

Total 

Nodes 

Gridl 

Surface 

Nodes 

Grid2 

Total 

Nodes 

Grid2 

Surface 

Nodes 

(«) 

114119 

1891 

51150 

1260 

(b) 

480G3 

1032 

21316 

516 

M 

C27G4 

1100 

25352 

550 


Table 1: Number of nodes used for (a) TAU0015, 
(b) TAU0015m, and (c) NACAOOlo airfoils. 

5 Two-Dimensional Results 

In this section, the fine-grid two-dimensional com- 
putational results are compared with experiments. 
A grid de-refinement analysis is then conducted to 
minimize the number of required grid points. The 
de-refinement analysis is necessary to minimize the 
cost of the three-dimensional computations. 

5.1 Fine Grid Comparison with Ex- 
periments 

All grids in the initial computations had the first, 
near- surface grid point below y + = 0.7 to ensure 
the sublayer of the turbulent shear flow was suffi- 
ciently resolved. Specification of the first grid point 
approximation is based on relationships among , 
Reynolds number, and the skin-friction for a flat 
plate boundary layer. For t urbul ent flow, this re- 
lationship is y + = A*/ * \/cj /2i?e, where cj ~ 
0.455/ In 2 (0'0GRe x ) and Ay is the physical distance 
for the first grid point away from the airfoil surface. 
Specifying a midchord y + = 0.102 will keep all near- 
wall points below 0.7 as will be shown later. 

The number of surface nodes and total grid nodes 
for the first comparison are shown in Table 1 for each 
airfoil. The TAU0015 airfoil grid has significantly 
more grid points because of the grid clustering near 
the surface discontinuities, which are not present in 
the TAU0015m and NACA0015 airfoils. 

Convergence of the steady- flow computations 
were achieved wffien the change in Cf w r as less than 
0.015% per iteration. At this criteria,, Cd changed 
by less than 0.005% per iteration. The only un- 
steadiness in the computations w’as observed at 
a = 26° where unsteady vortex shedding w T as ev- 
ident . 

Figure 4 show^s the computed C\ with variation 
in a for the XACA0015, TAU0015, and TAUOOISm 
airfoils compared with the experimental data. For 


the NACAOOlo, the maximum lift coefficient, C/. 
and stall angle, a, are 30%. and 4° higher than 
the experimental data. For the TAUOOlom airfoil, 
the maximum C/ and stall angle are 23%. and 2° 
higher than the experiments. Over predictions for 
the TAUOOISm airfoil are consistent with the earlier 
studies 16,17 which used the same airfoil but a struc- 
tured grid and two different RANS codes. In closer 
agreement, the computed results for the TAU0015 
airfoil approach the experimental results. The stall 
angle is overestimated by 1° and the maximum C/ 
by 15%. The subtle differences in geometry for 
the NACA0015, TAUOOISm and TAU0015 airfoils 
lead to large differences in the computed stall an- 
gle a and maximum C/. The drag coefficient com- 
parison indicates that all computations ovorpredict. 
drag prior to stall and underpredict drag in post- 
stall conditions. However, the TAU0015 results are 
in better agreement with the experiments than the 
NACA0015 and TAUOOISm results. 

A reasonable explanation for the differences be- 
tween the computational and experimental results 
for the TAU0015 geometry may be obtained with a 
careful review of the experimental results. In the 
experiments, (7/ and Cd are derived from numerical 
integration of the static pressure measurements at 
discrete locations, wffiereas the computations used 
pressure around the complete airfoil and included 
viscous drag contributions. Hence, a comparison 
of computed and measured pressure information is 
necessary for this review^. Figure 5 show’s computed 
pressure coefficients (C p ) for the TAU0015 airfoil 
compared with experimental results for a — 8° and 
a = 12° (pre-stall) and 22° (post-stall). The compu- 
tational results show’ pressure spikes at the leading 
edge and near 76% chord resulting from the geomet- 
rical discontinuities. Pressure spikes w’ere not mea- 
sured in the wind-tunnel experiments because taps 
could not be positioned on the actuator. Therefore, 
the experiments could not capture the additional 
pressure spikes predicted in the computations. 

Based on this understanding of the experimental 
data, Ci and Cd w T ere recomputed by integrating the 
pressure over the TAU0015 airfoil in regions consis- 
tent with the experimental pressure taps. That is, 
the contribution from the leading-edge actuator dis- 
continuity is now T excluded in the determination of 
Ci and Cd . The new T C( and variation w ith a 
are compared with experimental data in figure 6. 
The computed maximum Ci is now T overpredicted 
by 2% compared w ith the experimental data and the 
stall angles of attack are in agreement at 12% The 
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a 


Figure 4: Lift (Ci) and drag ( Cd ) coefficients ver- 
sus angle-of-attack (a) for experiments compared 
with computations for TAU0015, TAU0015m, and 
NACA0015 airfoils. 


overprediction in lift results in an underprediction 
of the drag. The post-stall (separated) conditions 
show notable disagreement between computational 
and experimental results. 

Because the SA turbulence model was not devel- 
oped or validated for separated flows, significant un- 
certainty exists in the computations for post-stall 
cases. Further, there is notable experimental uncer- 
tainty for post-stall conditions as confirmed in re- 
sults for a simple NACA0012 airfoil. 28,29 Hence, the 
computational and experimental results each have 
uncertainty in their respective quantities for highly 
separated flow conditions. 

5.2 Coarse Grid Solutions 

In any grid generation process, judgment must be 
made concerning the adequacy of the grids for the 
computations. This is usually accomplished by grid 
refinement studies. For the current study, the al- 
ternate approach of excessive grid refinement was 
used to obtain a solution that matches the experi- 
mental results. This approach allows assessment of 
confidence in the RANS methodology. To achieve 
the current grid for the TAU0015 airfoil, extensive 
grid clustering near the discontinuities was used to 
resolve any possible regions of flow re-circulation. 
In this section, the process of grid coarsening is dis- 
cussed. 

Table 2 shows the number of surface nodes and to- 
tal nodes for various grid coarsening strategies. The 
first row in Table 2 is the fine grid TAU0015 con- 
figuration shown as the first item in Table 1. The 
No Cluster case shows a decrease in the number of 
nodes by removing the grid clustering at the leading- 
edge, 76% chord, and the trailing edge. The No 
Cluster /Notch case combines removal of discontinu- 
ity grid clustering and removal of the notch geome- 
try at 76% chord. Coarse Grid 1 reduces the total 
number of surface nodes relative to the No Clus- 
ter /Notch case by approximately 2. Coarse Grid 2 
reduces the number of surface nodes on the actua- 
tor and lower surface of the airfoil by approximately 
2 compared with the No Cluster /Notch case. The 
finer grid on the upper surface is retained, however, 
so the separation point can be more accurately com- 
puted. 

Figure 7 shows lift and drag coefficients for the 
grid de-refinement cases compared with the origi- 
nal TAU0015 results. The coarse-grid results are 
sufficiently close to the TAU0015 experimental re- 
sults (within experimental uncertainty levels) and 
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Figure 6: Lift ( Ci ) and drag ( C<t ) coefficients ver 
sus angle-of-attack (a) for computations and exper 
iments for TAU0015 airfoil. 


Figure 5: Computational and experimental pressure 
coefficients ( C p ) for TAU0015 airfoil at a = 8° and 
a = 12° (pre-stall) and a = 22° (post-stall) for 
TAU0015 airfoil. 







Case 

Gridl 

Total 

Nodes 

Gridl 

Surface 

Nodes 

Grid2 

Total 

Nodes 

Grid2 

Surface 

Nodes 

(a) 

114119 

1891 

51150 

1260 

(b) 

69331 

1480 

30401 

712 

(c) 

68103 

1406 

30040 

716 

(d) 

36988 

716 

16439 

346 

(e) 

47965 

936 

21770 

482 


Table 2: Number of nodes used for (a) TAU0015, 
(b) No Cluster, and (c) No Cluster /Notch cases, (d) 
Coarse Grid 1. and (e) Coarse Grid 2 cases. 


Case 

Ci 

%Error 

(a) 

1.056 


(6) 

1.034 

2.1 

(c) 

1.030 

2.4 


Table 3: Lift coefficient (C/) for (a) Experiment, (b) 
TAU0015, and (c) Coarse Grid 2 cases at a = 12°. 
(Computational error relative to experiment) 

avoid the added expense of including the notch and 
grid clustering in the computations. Comparing the 
Coarse Grid 1 and Coarse Grid 2 results indicates 
that significant grid coarsening can be used on the 
lower surface, but additional coarsening on the up- 
per surface is not acceptable. Sufficient resolution 
on the upper surface is necessary to accurately com- 
pute the separation location. 

Results for the Coarse Grid 2 case are presented 
in a manner similar to that for the analysis leading 
to figure 6. The lift coefficient is recomputed for 
the Coarse Grid 2 case without the pressure con- 
tribution from the actuation region at the leading 
edge. Table 3 shows the TAU0015 and Coarse Grid 
2 solutions along^wltF the experimental result at 
a — 12°. (measurement station near maximum lift 
conditions). Over 65,000 points have been removed 
in the grid de-refinement process and the error in 
solution is approximately 2.4%. 

For the final series of two-dimensional computa- 
tions in this section, the effect of wall-normal grid 
resolution is examined to further minimize the re- 
quired number of nodes necessary for an accurate 
solution. In section 5.1, it was noted that in all 
previous grids, the first grid point was well below 
= 0.7 to ensure sufficient resolution of the sub- 
layer of the turbulent boundary layer. The spacings 




Figure 7: Lift (C/) and drag ( Cd ) with angle of at- 
tack for the de-refinement cases. 
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V 

y + 

Gridl 

Total 

Nodes 

Gridl 

Surface 

Nodes 

Grid2 

Total 

Nodes 

Grid2 

Surface 

Nodes 


0.102 

47965 

936 

21770 

482 


0.700 

36590 

912 

15573 

456 


1.382 

33455 

912 

13831 

456 


7.000 

27409 

912 

10442 

456 


Table 4: Number of nodes used for wall-normal grid 
coarsening. 

for the final series of simulations are shown in Table 
4. The resulting number of surface and total grid 
nodes are shown that result from increasing the dis- 
tance between the airfoil and the first wall-normal 
grid point. 

Figure 8 shows wall-normal grid coarsened lift 
and drag coefficients compared with the fine grid 
case. These results indicate additional grid coars- 
ening can be achieved in the wall normal direction; 
however, as expected, the viscous region can be- 
come under-resolved and the results begin to lose 
accuracy. Again, the initial spacing of the first 
wall-normal grid point is approximated using the 
flat-plate turbulent boundary layer relationship dis- 
cussed in section 5.1. However, a correlation for 
skin friction for airfoils does not exist to enable a 
determination of y+. For the cases in Table 4 that 
show relatively good quantitative agreement with 
the fine grid results, figure 9 shows the computed 
y + around the airfoil for a = 12°. The skin-friction 
varies with chordwise location (and a), therefore, 
y + varies with chordwise location. For these varia- 
tions in 2 / + , sufficient resolution overall is apparent 
with the good quantitative agreement in the results. 

5.3 Three-Dimensional Results 

In this section, three-dimensional pre-stall, stall, 
and post-stall flow conditions are investigated as 
a precursor to the active flow control simulations 
which will involve three dimensionality associated 
with actuator-induced flow fields. 

Ideally, the three-dimensional calculations should 
encompass viscous boundary conditions associated 
with the wind tunnel walls; however, such compu- 
tations are at this point prohibitly expensive. So 
initially inviscid walls with a subset of the wind 
tunnel domain are used for this first comparison 
of the three-dimensional computational results with 



a 


Figure 8: Lift (C/) and drag (Cd) coefficients ver- 
sus q for various near-wall initial grid spacing for 
TAU0015 airfoil. 



Figure 9: Computed with chordwise location for 
TAU0015 airfoil (a = 16°). 
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the two-dimensional simulation results and exper- 
imental data. Initial calculations were conducted 
in a volume grid that had a spanwise width equal 
to the airfoil chord. The outer boundary of the grid 
extended 10 chord lengths out from the airfoil geom- 
etry. Figure 3 shows the surface mesh at the lead- 
ing edge of the TAU0015 airfoil. This surface grid 
maintained the discontinuity at 76% chord. The 
minimum normal wall spacing on the airfoil was de- 
termined based on a =0.7, chord Reynolds num- 
ber of 1.2 million, and the skin friction for a flat, 
plate boundary layer. This grid had up to 26 layers 
of thin tetrahedral cells in the viscous region for the 
boundary layer with approximately 375 nodes along 
a chord wise surface cut. The total cell count for this 
volume grid was 1.85 million. Because the computa- 
tional results with this grid did not match the two- 
dimensional results or the experimental data, the 
grid was regenerated around TAU0015 airfoil with 
a much finer grid. 

As shown in figure 7 for the two-dimensional in- 
vestigation, grid resolution on the upper surface was 
necessary to capture the separation location and to 
accurately predict integrated aerodynamic perfor- 
mance. For the three-dimensional ’’fine” grid, the 
stretching in the viscous layer was modified to ob- 
tain up to 33 layers of thin tetrahedral cells for the 
boundary layer with approximately 800 nodes along 
a chordwise surface cut. With the increase in grid 
density, the span of the grid was reduced from 1.0 to 
0.5 chord lengths, resulting in a total cell count of 
2.9 million. The surface discontinuity at 76% chord 
was omitted for this three-dimensional grid because 
the two-dimensional results demonstrated that this 
discontinuity does not effect the aerodynamics of 
the wing and can be ignored for more efficient- com- 
putations. Figure 10 shows the surface mesh at the 
leading edge of the refined TAU0015 grid. 

Three-dimensional computational flow analyses 
were conducted on the TAU0015 for a Mach number 
of 0.15 and chord Reynolds number of 1.2 million. 
The angles of attack analyzed were S° } 10 n , 12° } 14°, 
and 22°. 

Convergence for the three-dimensional calcula- 
tions was determined by two methods. The first 
method involved tracking the solution residual until 
it dropped several orders in magnitude and started 
to level out. For the second method, the Q and 
Cd were monitored to verify that the solution had 
converged to a constant. The C p distributions were 
evaluated at the final solution and at 1000 iterations 
prior to the final solution to ensure they were not 



Figure 10: Three-dimensional fine grid case for 

TAU0015 airfoil. 

changing. Figure 11 shows the convergence history 
plots for the a = 8° case. The residual is reduced 
approximately 6 orders in magnitude and a near- 
constant lift coefficient is obtained in approximately 
2000 iterations. Convergence for these steady flow 
computations is considered attained when C{ and 
Cd changed less than 0.5% and 0.3%, respectively, 
over the final 1000 iterations. All the cases required 
between 3000 and 6000 iterations to achieve this 
convergence criteria, except for the a = 14° case 
which will be discussed later in this section. 

The variation of three-dimensional computed lift 
(C/) and drag (Cd) coefficients with angle of attack 
(a) for the TAU0015 airfoil are compared with the 
experimental data in figure 12. Both the coarse 
and fine grid three-dimensional computed results 
are shown. Here, the computed Ci and Cd are ob- 
tained by integrating the pressure distribution on 
the centerline of the model. The three-dimensional 
fine grid results are in better agreement with the ex- 
perimental data compared with the coarse grid re- 
sults, which are under-resolved. Futher discussion 
will involve only the fine grid three-dimensional re- 
sults. 

For low-a pre-stall conditions (a < 8°), the three- 
dimensional Ci and Cd results are slightly greater 
than the experimental results. At the higher-a, pre- 
stall region (a ~ 10°). the three-dimensional CJ and 
Cd results begin to more closely match the experi- 
mental results. Near stall angle conditions a = 12°, 
the computational results underestimates the max- 
imum Ci by 3% but matches the stall angle. For 
post-stall conditions (a = 14°), the computed C\ 
deviates from the experimental data; however, the 
computed Cd closely follows the trends of the ex- 
perimental data. With the separation bubbles on 
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Figure 11: Convergence of residual and lift coeffi- 
cient (Ci) with iteration for TAU0015 airfoil (a = 
8 °). 
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Figure 12: Lift ( Ci ) and drag ( Cd ) coefficients ver- 
sus angle-of-attack (a) for experiments compared 
with three-dimensional computations for TAU0015 
airfoil. 
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the upper surface becoming more developed at stall 
and fully separated flow at post-stall conditions, it 
is likely that the three-dimensional results are insuf- 
ficiently grid-resolved to capture the true flow char- 
acteristics and quantitatively match experiments. A 
second plausible reason for the near-stall and post- 
stall disagreement in the results resides with defi- 
ciencies in the turbulence model. The model was 
developed for attached and steady flows. 

Figure 13 shows the two- and three-dimensional 
computed pressure coefficients ( C p ) compared with 
the experimental data for a = 8°. Overall, very 
good agreement is achieved at this pre-stall angle 
of attack. Analysis of the three-dimensional com- 
putational results show the flow accelerating from 
the stagnation point around the leading edge re- 
gion of the airfoil and separating from the surface 
due to the backward facing step. The flow then 
reattaches at x/c= -0.002 followed by an acceler- 
ated flow region to x/c of approximately 0.01. The 
flow' then negotiates the increasing pressure gradi- 
ent at x/c=0.947 before separating at the trailing 
edge. Note, that the three-dimensional computa- 
tional pressure results overpredict both the tw'o- 
dimensional and experimental data near the leading 
edge region. 

Figure 14 represents a Mach contour map along 
the centerline of the three-dimensional grid for a = 
12° with the flow' separation region visible. The 
two and three-dimensional computed pressure coef- 
ficients are compared with the experimental data in 
figure 15. The three-dimensional results show' the 
flow' separating in the leading edge actuator region, 
with reattachment at x/c= -0.002 and a slight flow 
acceleration region existing back to x/c=0.006. The 
flow then negotiates the steep increasing pressure 
gradient back to x/c= 0.535 before flow' separation 
occurs back to the trailing edge. The calculated 
pressure rise from the three-dimensional results is 
just slightly steeper than the two-dimensional re- 
sults and the experimental data. Similar to the 
a = 8° case, the three-dimensional pressures are 
overpredicted compared with the tw'o-dimensional 
and experimental data. 

Considerable effort has been expended on the 
post-stall case of a = 14°. In order to capture to 
true flow' physics for this case, the original coarse 
surface mesh cell spacing w r as reduced along with 
increasing the number of cells in the boundary layer 
of the grid. Neither the coarse nor the fine grid so- 
lutions were successful in matching the experimen- 
tal aerodynamic performance of the TAU0015 for 



- 0.1 0.0 0.1 0.2 0.3 0,4 0.5 0,6 0.7 0.8 0-9 1.0 1.1 
x/c 



Figure 13: Tw'o- and three-dimensional computa- 
tional pressure coefficients (C p ) compared with ex- 
perimental results for TAU0015 airfoil at a = 8°. 
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Figure 14: Mach contour map at centerline cut for 
TAU0015 airfoil at a = 12° (insert of leading edge 
region). 


this angle of attack. The fine grid case was halted 
after 11650 iterations before the convergence crite- 
rion was met due to the slow but steady decrease in 
Ci. The C( at this point was 0.796. Examining the 
solution files showed the separation line progressing 
forward along the upper surface. Figure 16 shows 
the two- and three-dimensional computed pressure 
coefficient (C p ) for a = 14° compared with experi- 
ments. Clearly, the difference between the two- and 
three-dimensional results has increased for the post- 
stall conditions. However, in this comparison, the 
three-dimensional results are in much better agree- 
ment with the experimental data over most of the 
airfoil. 

The nearly converged surface pressure distribu- 
tion (after 11650 iterations) along with particle 
traces on the upper surface are shown in figure 17 for 
a = 14°. A distortion is observed in the separation 
front (spanwise curvature near leading edge). The 
particles traces from the leading edge run parallel to 
the flow field until the separation line is approached. 
A spanwise bias to the separation front is observed 
symmetrically on both sides of the airfoil centerline. 
In the reverse flow region, the particles traces move 
from just ahead of the trailing edge (x/c=.98) up to 
the separation line and again turn outward toward 
the walls. It may be that the in viscid boundary con- 
ditions are adversely impacting the flow such that 
the inviscid-side-w r all/leading-edge-juncture is caus- 
ing a spanwise bias to separation front. Similar re- 
sults were also found with the coarse grid solution. 
This behavior is not typical of inviscid boundary 




Figure 15: Two- and three-dimensional computa- 
tional pressure coefficients (C^) compared with ex- 
perimental results for TAU0015 airfoil at a = 12°. 
(lower graph is expansion of leading edge region) 
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conditions and remains to be resolved. 
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Figure 16: Two- and three-dimensional computa- 
tional pressure coefficients (C p ) for TAU0015 airfoil 
at a = 14°. (lower graph is expansion of leading 
edge region) 


Finally, figure 18 shows the two and three- 
dimensional computed pressure coefficients (C p ) 
compared with the experimental data for a — 22°. 
Reasonably good agreement is also achieved at this 
post-stall angle of attack over most of the airfoil. 
Three-dimensional computational results show the 
flow separating at the leading edge step with a small 
reattachment region from x/c= -0.001 to x/c=0.007, 
followed by a massive flow separation region back to 
the trailing edge. Both two- and three-dimensional 
result overpredict the pressure near the leading edge 
for this extremely stalled case. The flow separation 
can be visualized in figure 19. with the Mach contour 
map along the centerline of the surface grid. Results 



Figure 17: Three-dimensional surface pressure coef- 
ficient along with particle traces for TAU0015 airfoil 


at a = 14°. 


from the three-dimensional calculations, however, 
did not capture the small flow acceleration region 
back to x/c ~ 0.01 as found in the experiment . 

In the analysis presented in figure 6 for the two- 
dimensional comparison with experiments, remov- 
ing the contribution associated with the leading- 
edge actuator and re-integrating the pressures to 
obain Ci and Ca lead to good agreement between 
computational and experimental results. If the ac- 
tuator region is excluded from the integration in 
the three-dimensional analysis, the lift and drag de- 
crease (as expected). For a = 12°, the lift coefficient 
Ci decreases by only 1.1% while the drag coefficient 
decreases by 8.1%. As such, the computed drag co- 
efficient is in better agreement with the experimen- 
tal results when the actuator region is neglected. 
Also, the secondary leading-edge pressure spike over 
the actuator region in the three-dimensional compu- 
tations is smaller than in the two-dimensional com- 
putations; therefore, the change in C/ is smaller in 
the three-dimensional analysis compared with the 
two-dimensional analysis. 


6 Conclusions 

Results from two- and three-dimensional unstruc- 
tured grid Reynold- Averaged Navier-Stokes codes 
were compared against an experimental database for 
a NACA0015 airfoil modified at the leading edge to 
implement a fluidic actuator. This actuator caused 
a discontinuity at the leading edge. This paper doc- 
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Figure 18: Two- and three-dimensional computa- 
tional pressure coefficients ((7p) compared with ex- 
perimental results for TAU0015 airfoil at a = 22°. 
(lower graph is expansion of leading edge region) 



Figure 19: Mach contour map at centerline cut for 
TAU0015 airfoil at a = 22° (insert of leading edge 
region) . 


uments a comparison of computational results using 
the baseline airfoil (no control) with experimental 
data. 

The two-dimensional results show that the mid- 
chord discontinuity is insignificant and can be faired 
over for more efficient computational studies. The 
leading-edge discontinuity significantly affects the 
lift and drag; hence, the integrity of the leading-edge 
notch discontinuity must be maintained in the com- 
putations to achieve a good match with the exper- 
imental data. The pressure coefficients match the 
experimental data in all cases except for the highly 
stalled conditions. The comparison of computed 
lift and drag coefficients with angle-of-attack is in 
extremely good agreement with wind-tunnel data 
when the computational and experimental methods 
of computing the coefficients are the same. 

For all conditions, the three-dimensional pres- 
sures were slightly overpredicted in comparison with 
the two-dimensional and experimental data. The 
greatest overprediction occurred near the leading- 
edge of the airfoil. The three-dimensional, com- 
putational lift versus angle-of-attack results match 
the experiments in the pre-stall region (a ~ 10°) 
and the drag coefficient matches the experiments 
in post-stall conditions (a > 12°). For the near- 
stall and post-stall regions, the computational re- 
sults deviate from the experimental results, with a 
lower calculated C{ probably due to either lack of 
grid resolution or side-wall interference. The three- 
dimensional results show a spanwise bias to the sep- 
aration on both sides of the airfoil centerline; as yet 


15 



this bias is not completely understood. 

Future research will focus on refining the three- 
dimensional grid to more fully capture the baseline 
(no control) aerodynamics. The impact of far-field 
boundary locations will be assessed for separated 
flows relative to predicted performance. The cur- 
rent study will be extended to include RANS valida- 
tion of unsteady vortex shedding problems. Finally, 
actuator boundary condition models will be intro- 
duced to attempt the validation of unstructured- 
grid RANS for active flow control. 
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